Strong influence of the diffuse component on the lattice dynamics in Pb(Mgi/ 3 Nb 2 /3)03 
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The temperature and zone dependence of the lattice dynamics in Pb(Mg 1/ / 3 Nb 2 /3)03 is char- 
acterized using neutron inelastic scattering. A strong correlation between the diffuse and phonon 
scattering is measured. The lattice dynamics in Brillouin zones where the diffuse scattering is strong 
is observed to display qualitatively different behavior than that in zones where the diffuse scattering 
is weak. In the (220) and (200) zones, where there is a weak diffuse component, the dynamics 
are well described by coupled harmonic oscillators. Compared with SrTiOa, the coupling is weak 
and isotropic, resulting in only a small transfer of spectral weight from one mode to another. A 
comparison of the scattering in these zones to that in the (110) zone, where a strong diffuse com- 
ponent is present, reveals a strong coupling of the diffuse (or central) component to the acoustic 
mode. We speculate that the coupling to the central peak is the reason for several recent conflicting 
interpretations of the lattice dynamics based on data from zones with a strong diffuse component. 

PACS numbers: 77.80.-e, 61.10.Nz, 77.84.Dy 
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I. INTRODUCTION 

The relaxor ferroelectrics, with generic chemical for- 
mula PbB03, with the B site disordered, have generated 
intensive research recently due to their promising appli- 
cations as piezoelectric devicesi Despite much work on 
these materials there is no consensus on the origin of 
their interesting properties. Pb(Mg!/ 3 Nb2/3)03 (PMN) 
and Pb(Zn!/ 3 Nb2/3)03 (PZN) are prototypical relaxor 
systems which display a diffuse transition with a broad 
and frequency dependent dielectric response, but no well- 
defined structural transition without the application of 
an electric field. The nature of the diffuse transition 
and the low-temperature ground state is important to 
the complete understanding of these materials and has 
been the focus of much recent research and debate. 

Early studies of the refractive index by Burns and Do- 
cal (Ref. |2|) found that local regions of ferroelectric or- 
der were formed in a paraelectric background at a high- 
temperature, denoted as the Burns temperature T^. The 
connection of this high-temperature "transition" to the 
diffuse scattering (as observed by thermal neutrons) was 
made by Vakhrushcv et al. and Hirota et al. who ob- 
served that the diffuse scattering (around fru>=0 ) ap- 
peared at nearly the same temperature as Td~£ A cen- 
tral peak appearing at a temperature ~ T^ is one of 
the basic features of the random field model of Ref. 0. 
Recent work using cold-neutron elastic scattering (with 
much improved energy resolution) by Xu et al., Hiraka 
et al., and Gvasaliya et al. has characterized the diffuse 
component in detail&LS 

Pioneering work on the lattice dynamics was originally 
conducted by Naberezhnov et al. 9 Based on this, and fur- 
ther experiments around the (300) position, a new soft 
"quasi-optic" mode was proposed to explain the incon- 



sistency between the structure factors measured from the 
diffuse and the transverse-optic (TO) modern^ An alter- 
nate model was suggested by Gehring et al. and Waki- 
moto et aZiiiaiSii 3 . They suggested that the TO mode was 
the soft-mode and were able to describe the lattice dy- 
namics around the (300) and (200) positions consistently 
by a coupling between the TA and TO modes. Gehring 
et al. also observed a large dampening of the TO mode 
for q near the zone centre and attributed this to the for- 
mation of polar nanoregions. A different model was pro- 
posed by Hlinka et al. who were able to describe the 
lattice dynamics, including the anomalous dampening of 
the TO mode, in PMN qualitatively well using a stan- 
dard mode-coupling scheme, without the need for polar 
nanoregions^ 

The discrepancies in the interpretation of the lattice 
dynamics has been based solely on data taken in a (300) 
type zone, where a strong diffuse component is present. 
In all analyses of the phonons in PMN, there has been 
little consideration of the fact that the coupling may 
be anisotropic and also the feature that the (300) zone 
has a strong diffuse peak with an inelastic component. 
An example of anisotropic mode-coupling has been ob- 
served in SrTiOs by Yamada and Shirane M They found 
strong coupling between the acoustic and optic mode for 
phonons propagating along the [100] direction (measured 
near (200)), but little coupling was observed for phonons 
propagating along [110] (measured near (220))M> This is 
illustrated in Fig. which plots constant-Q scans taken 
near (200) in SrTiOs. At low-temperatures (4.5 K) a 
nearly complete transfer of spectral weight from the TA 
mode to the TO mode is observed, an indication of strong 
coupling. 

In this paper we describe neutron inelastic results 
aimed at elucidating the effects of mode-coupling (in- 
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FIG. 1: The upper panel schematically plots the effects of 
coupling as observed in SrTiOa for phonons along the [100] 
direction as measured in Ref. ITbI The lower panel displays 
the (HK0) scattering plane and illustrates the difference be- 
tween the Ti and T2 phonons. The dashed lines represent 
the regions in momentum space where the two modes were 
respectively studied. 



eluding any anisotropy) and the effect of the diffuse com- 
ponent on the phonons. Since all discrepancies to-date 
have been based on data obtained in Brillouin zones with 
a strong diffuse component (in particular (300)), we have 
focused our work in the (200) and (220) zones where there 
is no observable diffuse scattering. To study the effects of 
the diffuse scattering we contrast the results around (220) 
to preliminary data taken near the (110) zone, where a 
significant diffuse contribution has been measured around 



II. EXPERIMENT 

Neutron inelastic scattering experiments were con- 
ducted at the C5 spectrometer located in the NRU re- 
actor at Chalk River Laboratories. Constant-Q scans 
were made by fixing the final energy (E^-) to 14.6 meV 
and varying the incident energy (E^). A variable focusing 
graphite (002) monochromator and a flat graphite (002) 
analyzer were used with the horizontal collimation fixed 
at 12'-33'-S'-29'-72'. A graphite filter was used on the 
scattered side to filter out higher order neutrons and a 
liquid nitrogen cooled sapphire filter was used before the 
monochromator to reduce the fast neutron background. 
All inelastic data have been corrected for higher-order 
contamination of the incident beam monitor as described 
in detail elsewhere«i& 

The 9.3 cc crystal was grown by the modified Bridge- 
man technique previously described>iii& The sample was 
mounted in a cryo furnace such that reflections of the form 
(HK0) were in the scattering plane. The room tempera- 
ture lattice constant was measured to be a=4.04 A. To 
characterize further the sample we have measured the 
temperature dependence of the zone-centre transverse 
optic frequency and have found that it behaves identi- 
cally to that previously measured on much smaller sam- 
ples. As well, when plotted as a function of reduced tem- 
perature t = T/T c , the zone-center TO frequency agrees 
with that found Pb(Zn 1/3 Nb 2/ 3)0 3 £ 

We have focused our experiments in the (200) and 
(220) Brillouin zones where there is no observable dif- 
fuse component. The diffuse scattering (or sometimes 
referred to as the central mode) has been observed to 
be very weak in these two zones and therefore its influ- 
ence on the TA and TO modes is expected to be mini- 
mal. This approach differs from that in previous studies 
of the lattice dynamics which have addressed the prob- 
lem by comparing the dynamics in the (200) zone (where 
the diffuse component is weak) to that of the (300) zone 
(where the diffuse scattering is strong). 

As illustrated in Fig. we have focused our mea- 
surements on both the Ti and T2 phonons by conduct- 
ing constant-Q scans in the (200) and the (220) zones. 
In the (200) zone, we have studied transverse acoustic 
Ti phonons with the propagation vector along the [010] 
direction and polarization along [100]. For scans near 
(220), we have investigated T2 phonons with propagation 
vector along [110] and polarization along [110]. There- 
fore, any anisotropy in the lattice dynamics can be stud- 
ied by comparing these two zones without the influence 
of a strong diffuse component. 

To extract a linewidth, amplitude, and frequency we 
have fit the data to two uncoupled harmonic oscillators 
described by Lorentzian lineshapes as discussed in the 
Appendix (Eqn. The acoustic dispersion was approx- 
imated to be linear, w (q) = c|q|, where c is the phonon 
velocity. For optic modes, we have set the dispersion to 
have the form 
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FIG. 2: The dispersion of the TA and TO modes is plotted 
at 400 K and 100 K. The solid lines are the results of fits to 
the dispersion discussed in the text to extract the zone centre 
frequency of the TO mode. The 400 K dotted line indicates 
where the TO mode becomes overdamped and hence difficult 
to observe in a constant- Q scan. 
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FIG. 3: The TO zone centre frequency squared and the 
FWHM are plotted as a function of temperature. The open 
circles are from Ref. [l3| (around (200)) and the filled circles 
are data obtained in this study by extrapolating the measured 
dispersion relation to the (220) zone centre from higher q (as 
discussed in the text). 



^ (q) 2 = r»2 + a|q| 2 , (1) 

where f2 Q is the soft-mode frequency at q = and a is 
a constant for a given q directionii2iS£ For temperatures 
around T c we have used this dispersion to extract the 
zone center frequency tt by extrapolating from higher q, 
where the optic phonon is underdamped as discussed in 
the next section. 

For our coupled mode analysis we have used the same 
expressions presented in Ref. [2ll I22I Refs. [T2I and de- 
scribed in the Appendix (Eqn. [SJ. We have followed the 
coupling analyses of Refs. 0, |23, and Ref. by taking 
the coupling constant to be real. 

To extract the temperature dependence of the struc- 
ture factors we have used the formula for the intensity 
predicted by harmonic theory. The neutron cross-section 
as integrated over hui for a single phonon at w is ex- 
pressed as, 

^A|f(g)| 3 Q>Wll], (2) 

where | F(q) | is the structure factor at a particular mo- 
mentum transfer q and (n + 1) is the Bose factor. As 



measured in SrTiC>3, and as discussed in the next section, 
the temperature dependence of the structure factors is a 
sensitive measure of the coupling strength between the 
TA and TO modes. 



III. (220) AND (200) BRILLOUIN ZONES 

As stated above, to investigate the TA and TO 
phonons we have conducted constant-Q scans near the 
(220) and (200) positions (see Fig. 0J. The dispersion for 
both the TA and TO modes is summarized in Fig. |2Jand 
the temperature dependence of the full- width and the fre- 
quency position are plotted in Fig. Fig. 0] shows that 
at low temperatures the TO mode is underdamped while 
at higher temperatures the TO mode becomes heavily 
damped near the zone centre and is barely observable at 
400 K at g=0.1 r.l.u. Despite the fact that the TO mode 
becomes highly damped at q=Q.l, scans at q=0.15 and 
larger values of q show underdamped TO modes. This 
indicates that the anomalous damping of the TO mode 
around T c is strictly limited to values of q near the zone 
centre. 

To extract the zone centre frequency we have mea- 
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FIG. 4: Constant-Q scans at 100 K, 400 K, and 100 K, are 
plotted for q near (220). The solid lines are the results of fits 
to uncoupled harmonic oscillators. The parameters extracted 
from these fits were used to calculate the structure factors 
displayed in Fig. 5. 



sured the TO frequency position at large g, where an 
underdamped TO peak is observed, and extrapolated to 
lower q using the TO frequency dispersion described in 
the previous section. The results of the extrapolated zone 
centre frequency positions are plotted in Fig. [3] and are 
compared to those measured in Ref. [lousing constant-Q 
scans at the zone centre. The measured zone frequen- 
cies are also in good agreement with those obtained from 
infrared spectroscopy^ The agreement between the two 
sets of measurements is quite good and validates the ex- 
trapolation technique used in this work to extract the 
zone centre frequency. Concomitant with the softening 
of the TO zone centre frequency, we also observe an in- 
crease of the TO damping for modes near the zone centre. 
A similar broadening of the TO mode around T c has also 
been observed in PZN. 5 

A qualitative method of characterizing the effects of 
mode coupling is to study the temperature dependence 
of the integrated intensities as well as the frequency po- 
sitions. In the case of SrTi03 a large shift of spec- 
tral weight is observed to occur for phonons propagating 
along [100], but not [110] therefore indicating anisotropic 
mode coupling J£ In KTa03, strong coupling is also ob- 
served along the [100] direction as deduced by a large 
shift in the TA frequency with TO frequency^ Theories 



showing that that the observed frequency positions de- 
pends strongly on coupling have been described in Ref. 
Ill and Ref. IB 

Motivated by these previous studies we have charac- 
terized the effects of coupling by studying the phonon 
frequencies and the integrated intensities as functions of 
temperature. The temperature dependence of the fre- 
quency positions is illustrated in Fig. [2] where we show 
the dispersion of the TA and TO modes in the (220) zone 
at 400 K and 100 K. Despite a significant softening of the 
TO mode near the zone centre, the TA mode does not 
shift in frequency significantly. This result is confirmed 
by the constant-Q scans displayed in Fig. 01 which show 
little change in the acoustic frequency despite a large 
variation in the optic frequency. Based on the measured 
temperature dependences of the TA and TO frequencies, 
we infer that there are no strong coupling effects. 

To investigate the coupling near (220) and (200) more 
quantitatively, we have measured the structure factors 
from constant-Q scans for various values of q ranging 
from q=0.05 r.l.u to 0.3 r.l.u, examples of which are dis- 
played in Fig. To extract the peak positions, ampli- 
tudes, and linewidths of the TA and TO modes we have 
used the uncoupled harmonic oscillator lineshape con- 
volved with the resolution function as described in the 
Appendix. A nonconvoluted Gaussian centered at hu=0 
together with a constant term were used to describe the 
background. 

The temperature dependences of the structure factors 
provide a sensitive measure of the coupling. As the TO 
mode softens and the frequency gap between the TA and 
TO decreases for q near the zone centre we expect the 
effects of coupling to be dramatic. Specifically, as the 
two modes come closer together, a large transfer of spec- 
tral weight from one mode to another should occur in 
the presence of strong coupling. The amount of transfer 
of spectral weight will characterize the strength of the 
coupling. The structure factors are plotted in Fig. [SJas a 
function of temperature and compared with those mea- 
sured in SrTiOa. We have normalized all data by the 
structure factor measured when the TO mode is fully re- 
covered. In the case of PMN, the low-temperature struc- 
ture factor at 100 K was used while for SrTiOa the data 
were divided by the high-temperature value. This nor- 
malization puts the structure factors on the same scale 
and allows a direct comparison between different systems. 

The technique of characterizing the mode-coupling by 
tracking the integrated intensity has been used to in- 
vestigate the dynamics in several materials. The most 
relevant example is SrTiOs described in Ref. 15 which 
reports a significant anisotropy in the mode-coupling 
behavior based on structure factors^ For acoustic T2 
modes, propagating along [110], the integrated intensity 
was found to follow the prediction of harmonic theory. 
In the case of Ti modes, propagating along [100], the in- 
tegrated intensity was measured to deviate significantly 
(by more than a factor of 20) from harmonic theory (il- 
lustrated by filled circles in Fig. [SJ. Another exam- 
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FIG. 5: The measured acoustic structure factors, on a log- 
arithmic scale, as a function of temperature is plotted for 
PMN (open circles and solid line) and SrTiOa (filled circles 
and dashed line). SrTiC>3 shows a large change in the struc- 
ture factors for phonons propagating along [100], indicative 
of mode coupling. PMN, in comparison, shows only a mild 
change in integrated intensity. We have taken the harmonic 
structure factor to be the structure factor when the TO mode 
is fully recovered. 



pie where structure factors were used to characterize the 
coupling is the superconducting borocarbides described 
in Ref. |27] and Ref. |2^. These works report a large 
transfer of spectral weight between the acoustic and op- 
tic branches, emphasizing the large coupling between the 
two modes. The temperature dependence of the struc- 
ture factors in SrTiOs contrasts with that measured in 
PMN (open circles Fig. |5J); the latter shows a compar- 
atively small change in the measured structure factors 
with temperature (about a factor of two). 

For the Ti mode measured near (200) we find exactly 
the same behavior as that measured for the T2 mode. 
The (200) structure factors map directly onto the in- 
tensity measured around (220) when corrected for the 
Q 2 factor in the harmonic formula for the energy inte- 
grated intensity. We observe the same increase in spectral 
weight of the acoustic Ti mode with decreasing temper- 
ature as well as little change in the acoustic frequency 
position. We also observe similar symmetric lineshapes, 
well described by uncoupled modes, as presented else- 
wheieM^ Measurements similar to those described here 
have also been conducted in PZN near (220) and have 
found the same gradual increase in the acoustic inte- 
grated intensity below T c mtL In comparison to SrTiO;? 
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FIG. 6: Results of mode coupling fits to data at (2.0, -0.14, 
0) and (2.0, -0.21, 0) with the structure factors and coupling 
fixed as discussed in the text. The dashed lines are fits to 
uncoupled modes (A = 0), but with temperature independent 
structure factors 



and KTaC>3, the coupling between the acoustic and optic 
modes is very isotropic in PMN. 



IV. MODE COUPLING AND COMPARISON 
WITH THE (110) ZONE 

In PMN, the Ti and T2 acoustic modes are only weakly 
coupled to the optic modes as evidenced by the small 
change in the structure factors and the minimal shift in 
the acoustic frequency with temperature. This contrasts 
with the behavior in SrTiC>3 where the T2 mode shows 
little sign of coupling whereas the Ti mode shows evi- 
dence for significant coupling. The optic dispersion for 
both the Ti and T2 modes is also relatively isotropic. 
Measurements near the (200) position yield H 2 a=b.9 ± 
1.0 x 10 2 meV 2 -A 2 and near the (220) position h 2 a=4.6 
± 0.2 x 10 2 meV 2 -A 2 . Based on the symmetry of the 
dispersion and fact that the temperature dependence 
of the zone-centre TO mode measured near (200) and 
(220) agree exactly, we conclude that the TO mode is 
equally soft along the [100] and [110] directions. This 
also contrasts with the large anisotropy in both the dis- 
persion and temperature dependence observed in SrTiOa 
between the [110] and [100] directions. 
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FIG. 7: Constant-Q scans showing the temperature depen- 
dence of the T2 mode near the (110) position. An optic mode 
was not observable in our measurements indicating a small 
structure factor. The vertically dashed lines indicate the peak 
position of the acoustic mode measured in the (220) zone. 



We do observe an anisotropy in the linewidth of the 
transverse acoustic mode. For the acoustic Ti mode 
we observe a broadening of the acoustic mode above T c 
whereas for the T2 mode we observe no change in the 
linewidth around T c . We speculate that this anisotropy 
reflects the polarization direction of the correlated polar 
nanoregions at low-temperatures rather than anisotropy 
in the coupling. We emphasize that the broadening ob- 
served for the acoustic Ti mode is extremely small com- 
pared with that observed in the optic mode above T c ii2 
We are currently working to extend this result to higher 
temperatures. 

To characterize the mode-coupling in terms of the line- 
shape in more detail, we have fit the formula described in 
Ref. I21I and reproduced as Eqn. [SJin the Appendix. The 
coupled mode formula contains seven free parameters and 
therefore it is necessary to use other information to fix 
certain parameters. We have fixed the structure factors 
to those values which give a good fit over a broad range 
in g=0.1-0.25 r.l.u. and temperature T=100-500 K near 
the (220) and (200) zones. Based on this analysis we 

- 0.70 ± 0.15 



find that a value of | F opt i C /F t 



optic/ "acoustic 



describes the data reasonably well. This value is com- 
parable to that of 0.9, calculated based on the relative 
weight of Last and Slater modes discussed in Ref. 0. 
Along with q and temperature independent structure 



factors (F acoustic and F opt i C ) a temperature independent 
coupling constant was used which varied with q accord- 
ing to A = 2.7(±0.3) x 10 2 C 2 - 2.3(±0.5) x 10 3 C 4 meV 2 
(taking q=(C, C> 0)a*). The exact form of the coupling is 
defined and discussed further in the Appendix (see Eqn. 
|SJ. This form for the real coupling was previously used 
to describe the coupling in KTaOa^ We note that this 
q-dependent form for the coupling is exactlythe low q 
limit of that used in the calculations of Ref. [lj. 

Based on the profiles shown in Fig. 0] we do not ex- 
pect the coupling to be strong as the phonon peaks are 
well described by symmetric Lorentzians. Strong cou- 
pling has been previously observed to distort the phonon 
lineshapes dramatically due to phonon interference ef- 
fects. This effect has been observed in KNbC>3 in Ref. 
I21I and in BaTiC>3 in Ref. In both these examples, a 
symmetric Lorentzian lineshape gives a poor description 
of the phonon profiles. 

The results of these fits are illustrated in Fig. |S|for a se- 
ries of temperatures at (2, -0.14, 0) and (2, -0.21, 0). The 
agreement between the data and the model over a broad 
range of temperature and wave vector q is excellent. The 
dashed line in Fig. [fJlis the result of a fit to uncoupled 
oscillators with the structure factors fixed to the average 
value obtained from unconstrained fits at each tempera- 
ture (shown in Fig. [SJ . The agreement of the dashed line 
with the data is also very good showing only significant 
deviation at high temperatures. This analysis illustrates 
further that the effects of coupling between the acoustic 
and optic modes are relatively small. 

After fixing the structure factors and coupling we still 
find that a temperature and wave vector dependent Tto 
(with Tto increasing with decreasing q) gives a good de- 
scription of the data around T c . This implies that mode- 
coupling alone is inadequate to describe the anomalous 
damping of the TO mode in the "waterfall" region. In- 
stead, theories which involve coupling of the TO mode 
to some pseudospin (stochastic) variable 3 ^, presumably 
generated by random field o 5 i 31 , may be more appropriate. 
Such an extra coupling may capture the random nature 
of the polar nanoregions or any random fields introduced 
by chemical disorder. 

So far we have investigated the lattice dynamics in the 
(200) and (220) zones where the diffuse scattering has 
been shown to be weak or unobservable. To investigate 
the influence of the diffuse component we have conducted 
constant-Q scans in the (110) zone. The (110) Brillouin 
zone is the simplest region in which to investigate the ef- 
fects of the diffuse scattering as the optic structure factor 
is observed and calculated to be negligible. 

Fig. summarizes the results of the T2 phonon near 
the (110) position, a zone where a strong temperature 
dependent diffuse signal has been observed. In this zone 
we see a strong temperature dependence of the acoustic 
frequency position and lineshape, in contrast to that ob- 
served along the same direction in the (220) zone. It is 
difficult to reconcile this behavior in terms of coupling 
between the acoustic and optic modes as the change in 
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lincshape is much more dramatic than that measured in 
the (220) zone. Furthermore, the optic structure fac- 
tor has been measured to be extremely small. There- 
fore, its coupling should have no effect on the acoustic 
T2 mode measured near (110). Since the distortion of 
the TA mode in the (110) position cannot be attributed 
to the TO mode we speculate that there is a strong cou- 
pling between the diffuse component and the acoustic 
mode. This has recently been further illustrated by the 
cold neutron work described in Ref. Q. 

We speculate that it is this additional coupling to the 
central peak which is responsible for the anomalous lat- 
tice dynamics in Brillouin zones where a strong diffuse 
component is present. All interpretations of the lat- 
tice dynamics have been based on a coupling between 
an acoustic and optic mode and these conclusions need 
to be revaluated. We are currently working to extend 
our results to a (300) type zone where a strong optic and 
diffuse component exist. 

The importance of the diffuse component (or central) 
peak has recently been emphasized in Ref. in the con- 
text of interpreting relaxors in terms of order-disorder 
transitions. Our results, suggest that the diffuse com- 
ponent also plays a much more important role in the 
lattice dynamics than any coupling between the acous- 
tic and optic phonons as previously believed, and as is 
the case in conventional displacive ferroelectric systems 
like KTa03. It is interesting to note that the apparent 
softening of the acoustic mode near (110) occurs close to 
T c , and also in the same temperature region where the 
intensity and correlation length of the diffuse scattering 
have been observed to increased This is strongly sug- 
gestive that the coupling between the central component 
and the phonons is important in the relaxor transition. 

Coupling between a vibrational mode and a central re- 
laxational mode (diffuse component) has been observed 
in many systems and such a coupling can result in a soft 
mode analogous to a displacive transitioni^iiM For exam- 
ple, such a model has been applied to calcite which dis- 
plays both a soft mode and a continuum of scattering^ 
The strong correlation between the diffuse scattering and 
the acoustic mode in PMN may suggest that such models 
are suitable to explain the lattice dynamics behind the 
relaxor transition. They, of course, follow naturally from 
the random field models of Refs. and 0. We are cur- 
rently working to extend our low-energy acoustic phonon 
work to other zones to characterize this coupling and its 
symmetry in more detail. 



V. CONCLUSIONS 

The experiments presented here show that the lattice 
dynamics are very different between zones with a dif- 
fuse component and those where a central peak is ab- 
sent. This has been established by contrasting measure- 
ments made in the (220) and (200) zones with the (110) 
zone. The effects of mode-coupling have been charac- 



terized by studying the lattice dynamics in the (200) and 
(220) zones and have been shown to be weak and isotropic 
compared with those in SrTi03. We also find that mode 
coupling cannot exclusively explain the anomalous damp- 
ing of the TO mode observed in both the (200) and (220) 
zones. Further work is currently being pursued to study 
the dynamics in relation to the diffuse scattering in other 
zones. 
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VI. APPENDIX: PHONON CROSS SECTIONS 

In order to obtain detailed information on the 
linewidth from the inelastic spectrum, a model must be 
convolved with the resolution function and then fit to the 
observed scattering. In this section we state the cross- 
sections used to describe the constant-Q data presented 
in the main part of the text. 

Measurement of the neutron scattering intensity pro- 
vides a direct measure of S(q, u>), which is related to 
the imaginary part of the susceptibility x"(q,aj) by the 
fluctuation dissipation theorem^, 

S(q, W ) = 7r->H + l]x"(q,^), (3) 

where the energy transferred to the sample is defined 
by Tiuj=E % ~ E f and n(uj) = l/(e huJ / kBT - 1) is the Bosc 
factor. In the two sections contained in this Appendix we 
explicitly describe the two forms used for the uncoupled 
and coupled mode analysis. 



A. Uncoupled Harmonic Oscillators 

To extract the integrated intensity (and hence the 
structure factor) as a function of temperature and wave 
vector we have described the phonon cross section for 
a particular mode by the antisymmetrized linear combi- 
nation of two Lorentzians. The imaginary part of the 
susceptibility is given by, 
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C = wp?i(w2 - lo 2 ) + T 2 (u 2 1 - lo 2 )}, 



.4 



(4) 



[T(w) 2 + (hcu - ?kJ (q))1 

~ [r(w) 2 + (Tiw + ^o(q)) 2 ] ' 

where T(u>) is the frequency dependent half-width-at- 
half- maximum (HWHM), wo(q) is the undamped phonon 
frequency, and A is the amplitude. For acoustic 
modes, we have approximated the dispersion to be linear 
wo=c|q|. For optic modes, we have set the dispersion to 
have the form described in Eqn. [21 



B. Coupled-Mode Analysis 

To study the effects of mode coupling on the line- 
shape we have used the coupled mode formula previously 
derived and used to study perovskite ferroelectric sys- 
tems^i*^ For a pair of interacting modes, the imaginary 
part of the susceptibility can be written as 



X"(q,") 



(CY - BZ) 
(B 2 + C 2 ) ' 



(5) 



where, 



Y = F 2 (uj 2 



F 2 



lo 2 ) -2F 1 F 2 (lu 1 co 2 ) 1/2 A, 



B = (uj\- lo 2 )(lo 2 - lo 2 ) - w 2 rir 2 - WlW 2 A 2 



Z = lo(F 2 T 2 + F 2 T 1 ). 



The frequency positions and dampening constants for the 
two interacting modes are given by u>i, lo 2 , T\, and T 2 re- 
spectively. The structure factors for each mode are given 
by Fi and F 2 . The real coupling constant is defined by 
A. In our analysis we have followed previous works and 
assumed the coupling constant to be real. A generalized 
formula with both imaginary and real coupling constants 
can be found in Ref. |2J 

It is clear from the formula above that, in the pres- 
ence of coupling, a temperature dependent frequency will 
result in an apparent temperature dependence (as mea- 
sured from the integrated intensity of the two phonon 
peaks in the intensity profile) in the structure factors. 
This was illustrated in the early work of Ref. and 
provides the basis for our fits over a broad range in tem- 
perature and energy transfer to the phonons in the (200) 
and (220) zones. 
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